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ABSTRACT 

Using standard big-bang nucleosynthesis and present, high-precision measurements of light element abun- 
dances, we place constraints on the self-gravity of radiation pressure in the early universe. The self-gravity 
of pressure is strictly non-Newtonian, and thus the constraints we set are a direct test of this aspect of general 
relativity. 

Subject headings: general relativity; Friedmann equation; big bang nucleosynthesis 

1. INTRODUCTION 

Certain aspects of general relativity are well tested. For ex ample, the Schwarzschild metric h as been quantitatively verified 
in the weak-field limit on small scales, e .g., the Solar system ( Shapiro|[T964l ; iBertotti et alj |2003') and binary radio pulsars (e.g.. 



Hulse& Taylor 19751 iTaylor et alJfT979 : Weisberg & Tayloiill984 r and on galaxy scales (e.g., Bolton, Rappaport, & Burle^ 



20061). In another fundam ental test of g eneral relativity, the existence of gravity waves has been established (e.g., 
jTavlor & Weisb erg 1982; Weisberg & Ta vlor 1984, 2003). Gen eral relativity theory, in the form of the Fried mann-Robertson- 
r-| I Walker rnetric (lFriedmann| ll922L 1924: Robertson 1935, 1936: IWalkeHl936h and the Friedmann equations (iFriedmannll 19221: 
|. lLemaitrdll927h which govern the expansion behavior of the Universe are used extensively in cosmology, and are at the core of 
this carefully and tightly woven paradigm. However, it is probably fair to say that the Friedmann equations, while providing a 
very self-consistent and highly successful framework for cosmology, have not been subjected to extensive, independent testing. 
In this paper we show that one particular aspect of general relativity, i.e., the self-gravity of pressur e can be tested quantitatively. 

The development of sophisticated bi g-bang nucleosynthesis (B BN) codes dWagoner et al.llT967b . coupled with measurements 
of the relevant nuclear reaction rates (iCaughlan & Fowled 1 19881) , have allowed observations of light elemental abundances to 
become powerful tools with which to investigate the evolution of the early universe. Computational predictions over a wide 
range of parameter space, when compared with observati ons, have yie lded constraints on the curr ent-epoch baryon density 
^ ! JWa g oner 1973) , neutr ino physics (Yang et al. 1979; KawasakTet al.l 19941) , the fine str ucture constant (Bergstrom et al. 199 9), the 
■ gravitational constant (lYang et aLlll979l : ICopi et al.ll2004l) . primordial magnetic fields (iKernan et al.lll996l) . and other parameters 
Jyr ' of astrophysical interest. 

Increasingly exact measurements of elemental abundances, as well as augmented understanding of the processes (i.e., stel- 
lar nucleosynthesis) which hav e altered the original abundances, allow these restrictions to be continually refined. Deuterium 
abundances (Black & Dalgarno 1973; O'Meara et al. 2006) , helium abundan ces (Peebles 1966; Izotov et al. 2007), and lithium 
abundances, including the effects of stellar mixing ( Pinso nneault et al]l2002h have all been well measured. More recently, ob- 
• servations of the co smic microwave background (CMB) have yielded an independent estimate of ?/, the baryon to photon ratio 
V$ ■ (ISpergeletalJl2007h . 



2. ANALYSIS 
2.1. Friedmann Equations 



H 

' In the standard Friedmann-Robertson-WaUcer (FRW) cosmology, the expansion of the scale factor, a, of the universe as a 
function of its density and pressure is given by the Friedmann equations. In terms of a, the "acceleration" of the scale factor is 
given by: 

a AttG ( 3P' 



P+— ■ (1) 

where pc^ is the energy density and P is the pressure. This is an exact solution of the Qrr Einstein field equation for a homogenous 
and isotropic universe. Note that the 3P term, implying the self-gravity of pressure, has no analog in Newtonian gravity. The 
other well-used form of the Friedmann equations: 

—P , (2) 

arises from the 0„ component of the Einstein field equations (for a flat Universe). For any fluid with a significant pressure, and 
given its appropriate equation of state, eq. (1) can be integrated to yield eq. (2) - but only, of course, if the 3P term is included. 

It could be argued that eq. (2) is adequate to derive the evolution of the scale factor, and yet it depends only on p, and not 
explicitly on P. Yet, we know that without the 3P term in eq. (1), the constant of proportionality in front of p in eq. (2) will be 
incorrect, as long as pressure is significant. The Friedmann equations have been frequently derived fr om Newtonian theory, in 
the context of a pressureless (P = 0) universe, both for scientific and strictly pedagogical reasons (e.g., iMilne & McC"realfT93l 
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Fig. 1. — Combination of tlie various constraints on x- Tlie red-liatched region indicates tlie independent WMAP constraint on rjio (i.e., the baryon/plioton 
ratio in units of lO'"). The bold contours represent the observed estimates and their 1 a errors. Green contours are deuterium abundances (D/H)p. Blue contours 
are for the mass ratio of He-4 (Yp). Purple contours represent lithium abundances (jLi/H)p. The black ellipse is the 90-percent confidence contour based on the 
helium-4 and deuterium measurements. 

lUzan & Leh oucq 2001; see also lJordanl2005l) . However, since the concept of self-gravity of pressure does not exist in Newtonian 
gravity, it follows that eqs. (1) and (2) cannot be properly derived from any purely Newtonian construction CRainsford 2000,) . 
Finally we note that the Friedmann equation given by eq. (2) (in which P does not explicitly appear) formally "knows" about the 
form of eq. (1) (which does contain P) via the Bianchi identity. 

Consider now the radiation dominated epoch of the universe which prevailed for the first few thousand years after the big bang. 
Let us start with the Friedmann equation in the form of eq. (1). From special relativistic considerations alone, we know that 
P=\ /3>pc^. Therefore the a equation can be written as: 

d 47rG(l + x) 

-= 5 P , (3) 

a 3 

where x is used as a placemarker for the pressure term (i.e., x = 1 if the pressure term is present, and x = if not). Since p oc 
for radiation and relativistic particles (e.g., neutrinos) undergoing an adiabatic expansion, eq. (3) can be integrated to yield: 

SttG/I + X 



(4) 

Variation of the value of x away from unity would would significantly affect the outcome of BBN, degenerate with a variation in 
G (the gravitational constant) or Ni, (the number of neutrino families). However, assuming these latter two to be constant at their 
nominal values of 1 and 3 respectively, nucleosynthetic constraints are able to place limits on x and thus test the validity of this 
aspect of general relativity in this regime. 



2.2. Nucleosynthesis Calculations 

Nucleosynthesis calculations were performed with the canonical Kawano BBN codes (iKawanoll 1992b which have been modi- 
fied and updated with the latest reaction rates. The codes offer the ability to alter relevant early universe parameters such as rj and 
G. Since the parameter we seek to constrain, i.e., (1 + x)/2, is multiplicative with G, we simply vary G as a surrogate. Altering 
the parameters over repeated runs of these codes yields abundances of light elements as a function of rj and x- 

In this work ^75000 BBN calculations were performed, varying x between -1/2 and 7 in steps of 0.02 (or equivalently G/Go 
from 1/4 to 4 in linear steps of 0.01) and varying logjg -q from -10 to -9 in steps of 1/200 dex. The result of these calculations, is 
a grid of {rj, x} pairs, each point with a set of associated elemental abundances (see Fig. 1). 

As estimated in B uries e t al. (2001b) the uncertainties in the nuclear reaction rates leading to deuterium are 3.4% (95% confi- 
dence limit), and as estimated in lBurles et alj ( 1200 lal) >0.1% (95% confidence limit) for He-4. 
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2.3. Light Element Observations 

Deuterium provides an excellent constraint because, as has been understood for the last fe w decades, the obse rved amounts 
necessitate that it must have formed in the big bang rather than in stellar or ga lactic processes dReeves et al.lll97 J). It is equally 
fortunate that it is extremely sensitive to the baryon to photon ratio, rj ( Buries et al.l 1200 lab . Deuterium measurements from 
along the li ne of sight to high red shift quasars have culminated in the current determination of logjy(D///)p = -4.55 ± 0.04 (la 
confidence: IO'Meara et al.lF200 6'). Contours of constant deuterium abundance are shown as green curves in Fig. 1. 

Lithium-7 was initially a frontrunner as a nucleosynthetic constraint. However it h as somewhat fallen out of favor because 
of large uncertainties arising from the systematic effects in stellar depletion models. iPinsonneault et al] (12002!) determine the 
primordial lithium abundance as logj,)(L////)^ = -9.6 ± 0.2 (la confidence), which while consistent at the 1-a level with the 
estimates derived from other light elements, is unable to provide any additional significant constraint due to the relatively large 
error bars. Until this situation can be improved, lithium-7 will not likely play a strong role in constructing BBN-based parameter 
limits. Contours of constant lithium abundance are shown as purple curves in Fig. 1. 

Helium-4, though not as sensitive to rj as deuterium, also serves as a useful constraint on x- The primordial abundance is 
determined from present-day observations of HII regions, coupled with models of gal actic chemical evo lution. The primordial 
He-4 mass fraction, Yp, is most recently reported to be 0.2472 ±0.0012 (1 a confidence: llzotov et al.ll2007h . Contours of constant 
helium-4 abundance are shown as blue curves in Fig. 1. 

2.4. WMAP Constraints 

Observations of the cosmic microwave backgroun d are able to provide an independent estimate of the baryon/photon ratio r/io 
(i.e., 77 in units of 10'"). The WMAP collaboration (ISpergel et al.ll2007h gives as the three-year measurement rjio = 6.116;I;o 249- 
This measurement assumes x = 1' therefore does not provide any additional constraint. However, as shown in Figure 1, the 
WMAP value of 7710 is in excellent agreement with the light element constraints. 

2.5. Constraints on x 

Figure 1 displays the results of our analysis. A section of the x} parameter space is shown, with number density (relative 
to hydrogen) contours for deuterium shown in green, and for mass fraction of He-4 shown in blue. The purple contours represent 
the Li-7 number density (relative to hydrogen). The boldest contours indicate the current best abundance determinations, with 
the adjacent intermediate-width contours indicating the corresponding ranges of uncertainties (la confidence). 

Assuming statisically-independent Gaussian errors, one can calculate the probability, via a maximum likelihood analysis, 
that the abundance determinations agree with the corresponding results of the BBN calculations at a given point in the {77, x} 
parameter space. Only the constraints due to D and He-4 were used to compute the uncertainties in x and rj (see §2.3 for an 
explanation of why Li-7 was not included). The black ellipse is the 90% probability contour 

3. CONCLUSIONS 

As illustrated in Figure 1, the combined constraints correspond well with the general relativity prediction of x = 1 and the 
WMAP constraint on 7710. The 90-percent confidence contour gives a limit of x = 0.97 ±0.06. Thus, we have used current light 
element observations and BBN computations as a test of the general relativistic self-gravity of pressure. The agreement is good 
to within ~6%. We reiterate that our result can also be expressed as a constraint on the Newtonian constant, G (during the BBN 
epoch), the number of light neutrino families, or degeneracies among neutrinos (see §1 for references). Assuming that these 
latter parameters take on their nominally expected values, then the self-gravity associated with the pressure of light and neutrinos 
during the BBN epoch has been quantitatively measured. 

JS acknowledges support from the Paul E. Gray (1954) Endowed Fund for UROP. We thank Al Levine, Scott Hughes, and Ed 
Bertschinger for very helpful discussions. 
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